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Abstract. We investigate FRW cosmological solutions in the theory of modulus field coupled to 
gravity through a Gauss-Bonnet term. The explicit analytical forms of nonsingular asymptotics 
are presented for power-law and exponentially steep modulus coupling functions. We study 
the influence of modulus field potential on these asymptotical regimes and find some forms of 
the potential which do not destroy the nonsingular behavior. In particular, we obtain that 
exponentially steep coupling functions arising from the string theory do not allow nonsingular 
past asymptotic unless modulus field potential tends to zero for modulus field -0 — > ±oo. Finally, 
the modification of the chaotic dynamics in the closed FRW universe due to presence of the 



Gauss-Bonnet term is discussed. 
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A great progress in constructing a theory of all physical interactions occurred over the 
past few decades. Standard model of the particle physics unified electro-magnetic, weak 
and strong interactions into one theory at the energy scale of 10 16 GeV. Unfortunately, 
methods of Quantum Field Theory do not result with a huge progress when applied 
to the quantization of General Relativity. Therefore to avoid problems in cosmology 
and black hole physics a new approach (String Theory) was constructed to study the 
topological structure of the space-time in the near Planckian regions |2|. 

In the perturbational approach developed at a first stage string theory predicts 
the Einstein equations to be modified by the higher order curvature corrections in the 
range where curvature of space-time has near-Planckian values (the situation which took 
place at very early stages of our Universe evolution). At the present time the form of 
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these higher order curvature corrections in the string effective action is not investigated 
completely. We do not know the general structure of the expansion and, hence, the 
direct summation is impossible. But in the perturbative expansion the most important 
correction is the second order curvature one, it is the product of the Gauss-Bonnet 
with dilatonic and moduli terms. In four dimensional space-time Gauss-Bonnet term 
represents a total divergence but its combination with a dilatonic (moduli) field makes 
the contribution of the second order correction to be dynamical. The investigations in 
the frames of the discussed model were performed formerly and a set of new cosmological 
and black hole type solutions was found [§, ||, ||, § which provide new types of space- 
time topology 0]. All these effects appear in the higher order curvature string gravity 
and are absent in the minimal Einstein gravity. 

At the next stage of the string theories evolution it was proved that all five 
independent string theories were connected by the different dual transitions and they 
all were a part of M-theory || , which in low energy limit gives rise to 1 1th dimensional 
supergravity. Using duality transitions one has an opportunity to work with a non- 
perturbative theory but at the current moment non-perturbative M-theory is only at 
the beginning of its way, its study methods are not completely created |§. 

So, in order to make a little step in understanding the space-time structure, it 
is possible to work in the perturbational approach and try to study new topological 
effects of low energy string action with the higher order curvature corrections. Here it 
is necessary to note that all the conclusions obtained in the near Planckian region must 
be treated as only preliminary directions on some effects, later they have to be verified 
by quantum gravity calculations. So, the most convenient form of the string gravity 
action is (we use the units where m PL /^/8ii = 1): 



S =2 



(1) 



where R is Ricci scalar, -0 is modulus field, (f) is dilaton, A (<5) is dilatonic (modulus) string 
coupling constant. Their values depend upon the concrete type of the string theory. 
Second order curvature correction (Xe^+d^ScB represents the product of the coupling 
functions e~ 2< ^ and £(?/0 and Gauss-Bonnet combination Sgb — RijkiR l ^ kl — ^RijR 1 -* + R 2 ■ 
The importance of moduli field in cosmology is a consequence of the fact that the sing 
of S is not fixed from the first principles and can be either positive or negative depending 
upon the detail structure of the theory. Antoniadis, Rizos and Tamvakis |10[ have 
shown that in the case of negative 5 a new class of cosmological solutions appeared with 
nonsingular past asymptotic and smooth transition to ordinary Friedmann Universe. 
In the next papers 0, [TTJ, |T^] such kind of solutions was constructed for metrics with 



zero or positive spatial curvature. We recall that in Einstein theory with a scalar field 
all the possible past asymptotics for flat models are singular (thought if scalar field 
potential is very steep, trajectory can approach a singularity in a complicated way, see 
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I6|). In closed models there are nonsingular periodical and aperiodical solutions but 



their measure in the initial condition space is equal to zero [17]]. The presence of modulus 



term in action (|1]) leads to appearance of a set of nonsingular solutions with finite (and 
sufficiently large) measure of required initial conditions. 

The presence of dilatonic field does not change qualitatively these non-singular 
solutions. The authors of Ref. |Tl[ have proposed to neglect the dilatonic contribution 



at all to simplify the analysis and, hence, to show the phenomenon of nonsingular 
cosmological behavior induced by modulus field in its "pure" form. We follow this 
proposal and consider only modulus-dependent part in addition to the Einstein term in 
(HI). £(V0 m string theory is defined in terms of Dedekind function and is exponentially 
steep for ip — > ±00 ||. We do not restrict ourself by this particular case and analyze a 
much wider class of coupling functions. 

The concrete form of nonsingular asymptotic depends on the function £(V0- 111 [0 
the particular case £ ~ if) 2 is investigated both by the analytical and numerical methods. 
The asymptote for the closed Friedmann-Robertson- Walker (FRW) metric has the form 




const = V-L25 , 7 p = < l p + J- t, (2) 



where a is the scale factor of the FRW metric. 

One can see that the nonsingular dynamics on a is accompanied with the growth of 
the absolute value of the modulus field. The same feature presents in the FRW flat case 



see below). It was noticed in [fn}| , that in such situation the modulus field potential, 



omitted in ([]]) can be important. As nonzero modulus field potential appears in theories 



with supersymmetry breaking |14| and the scale of SUSY breaking is very small in the 
Planckian scale, therefore, modulus field is very light. But the unbounded growth of 
modulus field causes the necessity to take into account this potential term in asymptotes 
(2), and such puzzle requires special additional analysis. So, in this paper we consider 
the action in the form: 



S = — / <l l .\-\f^7j 



-R + d^d^ + 5^)S GB + 2 V{$) 



(3) 



where V(ip) is the modulus field potential. 

The structure of our paper is the following: in Sec. 2 we present nonsingular 
asymptotics in FRW flat case for several types of modulus coupling function 
and investigate the influence of the modulus field potential upon them. A particular 
attention is paid to possibility of the nonsingular regime to survive in the presence of 
the nonzero potential. In Sec. 3 this analysis is extended in a more qualitative way to 
the positive curvature FRW metric. Sec. 4 provides a summary of the results obtained. 
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2. Flat case 

In this section we consider dynamics of a flat FRW universe with the metric 

ds 2 = dt 2 - a 2 [t) [dx 2 + dy 2 + dz 2 ^j . (4) 
After introducing a Hubble parameter h = a/a the field equations are 

3h 2 (l + A5^ -ij 2 /2-V = 0, (5) 
2(h+ h 2 ^j (l + 45£h) + h 2 (l + 46£) + 4> 2 /2 -V = 0, (6) 
i) + 3hip + V - 125£'h 2 (h + h 2 ) =0 . (7) 



The prime (dot) denotes differentiation with respect to if) (t). 

In this section our aim is to construct the explicit form of nonsingular asymptotic 
solutions of the equations ©-(0) • In order to attain this aim we consider the contraction 
of the model, therefore, the expansion of the flat universe corresponds to the time reverse 
of the solutions listed below. 



2.1. The case V{ip) = 



Previously Antoniadis, Rizos and Tamvakis [|10| , |TTJ have shown that nonsingular 
solutions can exist only for 5 < 0, which would be assumed from now on. 

In our investigation we assume a modulus coupling function to be represented by 
a power law — ip n+1 /( n + 1)) or asymptotically exponential ~ & ) forms. 

We do not consider functions £(V0 growing slower than if) 2 because for — > oo these 



functions do not provide the violation of both strong and weak energy conditions [[!(], [T2|] . 

To construct nonsingular asymptotics it is possible to neglect the first order 
curvature term in ([5]) in comparison with the second one. So, the reduced constraint 
equation is 

^ = 2A5h 3 . 

The left hand side (LHS) of the reduced constraint equation can be easily integrated 
for considered functions £(V ; )- The result is 

LHS ~ { , C(^) ~ ^ n , n>2; 




Thus, these cases must be treated separately. 



For £(if>) = V' 2 /2 anticipating the asymptotic solution 
h = At a , ^ = Be pt 
and substituting this ansatz into the reduced field equations one obtains 



h(t) 



V485 ' 

ip{t) ~ A exp (24<5/i 3 t) , A = const. (8) 

The scale factor a(t) in this case exponentially tends to zero. 

For £(ip) = ?/>( n+1 ) /{n + 1) anticipation of the asymptotical solution in a form 

h = At a , ip = Bt 13 

leads to two different cases: 1 < n < 5 and n > 5. For the first one the explicit form of 
the solution is: 

, / x n — 5 i 

hit) ~ ; r* , 

w 5(n-l) 

/l25(n- 1) 2 \^ 2 

re — 5 

the scale factor on this solution tends to zero vie power law a(t) ~ t . For the 
second one we have: 

2n 

~ At ™+s , A = const < , 

i 

(5 \ u_1 5 

24^ + 5)^ ) (10) 
For n > 5 the scale factor a(t) asymptotically tends to a constant positive value. 

The case £,(ip) ~ e'^', which is the asymptotic form of function £,(ip) arising in 
string theory has been previously investigated in detail in |10| . In the cited paper three 
possible asymptotical solutions have been listed but only one of them corresponded to 
the contraction of the Universe. In order to make a complete reference of nonsingular 
asymptotics we specify it in our notations here: 

h{t) ~ At~ 2 , A = const < , 
\m\ ~ 51nt + ln^^. (11) 

2.2. The case V{^) ^ 

Here it worth to notice that for arbitrary functions £,{ip) and V(ip) the action ([3]) allows 
rescaling with respect to 5: 

i — h V 

t -> M\t , /!_>", V -> -= . (12) 
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Hence, it is possible to eliminate the parameter S from the equations by 
substituting S = — 1 which will be assumed from now on. 

First of all it is necessary to check whether it is possible to introduce the potential 
V(if>) which does not violate the asymptotical solutions obtained in the previous 
subsection. Substituting these asymptotical solutions into constraint equation (^) one 
can see that all the items in constraint tend to zero with t — > oo unless £(V0 = 
if) n+1 /(n + 1) for n < 3. 

For £(ip) = %- the asymptotical solution in the form (Q) provides all the terms in 
(I) to be of second power with respect to if). So, the nonsingular solutions are not 
violated if the potential grows slowly than if 2 at \if>\ — > oo. For the particular potential 
V(ip) = m 2 ip 2 /2 a generalization of the solution (||) (h = const , if) = Ce At , A = const , 
C = const) can be substituted into the field equations which leads to 

A 2 + 24h 3 A + m 2 = , 
A 2 + 3hA + 12h 4 + m 2 = . 

This system has real solutions only for m < m w 0.32 which implies that the 
nonsingular asymptotes are violated unless m < mo. The potentials steeper than if) 2 
always violate the nonsingular asymptotic under consideration. 

For £(if)) = the asymptotical solution (Jfty gives if) = const. So, the only 
possible type of potential which does not violate this nonsingular asymptotic is the 
asymptotically flat one {V(if>) —>■ A = const for \if>\ — > oo). Substituting the ansatz 
h(t) ~ Ct , if)(t) ~ At (C = const, A = const) into the field equations the constraint 
gives C 3 = — • Being combined with the equation (|7p it results 

(2 A - 5 A 2 ) 3 + 24A 4 (A 2 + 2 A) 2 = , 

which has real roots only if A < Ao ~ 0.42. This means that the solutions remains 
nonsingular for A < Ao- 

As it is stated above, for £(if>) = ifs n+1 /(n+1), n > 3 and for exponentially steep £(ip) 
all terms in (|5|) tend to zero which implies that an arbitrary asymptotically nonvanishing 
modulus potential destroys corresponding solutions flID"D,(|TID. 

At the end of this section we describe singular solutions of our system. Previously it 
was shown [IT] that for V(if>) = there was only one singular asymptotic. It belongs to 



the class \h\ — > 00, if) — finite of singular solutions. It turns out that the introduction of 
a positive potential V(if>) which remains finite with it's first derivative for finite if) will 
not lead to new singular asymptotics of this type. To prove this we apply the procedure 
described in the above cited paper. 

The equations of motion for the system (0)-© may be rewritten as: 



if} = l25£h s ±JU45 2 £> 2 h 6 + Gh 2 -2V, (13) 



dip ipB 
where 

A = 48^> 4 /i 4 + 5ip 4 h 2 + 36h e + 36ip 2 h 4 + h 2 V 2 + 12h 3 ipV - 2ip 3 hV 

-2V (l2/i 4 + Qip 2 h 2 + 2hipV') , (15) 

B = 3Qh 4 - Y2h 2 ip 2 + 5^ 4 + 12V (ip - 2h 2 ) + AV 2 . (16) 

The equation (|7|) can be excluded since the equations are not independent. 

For the singular solution under consideration, keeping in mind finiteness of ip and 
aforementioned restrictions on V(ip), we can neglect the term V(ip) in equation ([l"3|) 



and, therefore, we only need to study equation (|14|). There are three possibilities: 
5£'h 2 ~ 0(1), 5£'h 2 < 0(1) and 5£'h 2 > 0(1). In the first case we obtain from (0) 
that ip ~ ah, a = const. Substituting this solution into equation fll5|) we can see that it 
is not necessary to proceed the separate cases for It is enough to notice that we can 
neglect all l/-terms (the eldest of them are proportional to h 4 ) with respect to terms 
proportional h 6 . In equation (|l(J) we also can neglect all \^-terms (the eldest of them 
are proportional to h 2 ) with respect to terms proportional to h 4 . Thus, this case can 
be treated as with V = 0. In the second case we obtain ip ~ ±y/Eh and we treat it in a 
complete analogy with the previous one. For 5^'h 2 ^> 0(1) there are two possibilities: 



ip + ~ £'/i 3 and ip- ~ Jr. For ip = ip + the eldest l^-term in equation flT5]) is proportional 



to h 10 which can be neglected with respect to h 1A terms. For ip = ip_ the eldest V-teim. 
in ([TB]) is proportional to h A which can be neglected with respect to h 6 term. 
So, the asymptotic ]ll 



h(t) ~ (t — t ) 1 — > — oo , t — ► t — , t = const , 

1 2 

W) ~ 4>o - (t - to) , 4>o = const (17) 

does not depend on V(ip) if V(ip) satisfies the aforementioned requirements. 

If V — the signs of h and ip are crucial as they are conserved independently. For 
V 7^ the sign of ip is not conserved. So, solution can be singular irrespectively of the 
initial sign of ip in contrast to the case V = 0. 

3. The positive curvature case 

The FRW metric for the positive spatial curvature case is 

ds 2 = dt 2 -a 2 (t)dtt 2 , 
where dfl 2 is the volume element of the 3-dimensional sphere, field equations are 

2-(l + + (~ 2 + ±) (1 + 4$ = V(iP) - % , (18) 



-3^- 125^ (~ 2 + 4) + = , (19) 



a a \ cr af 



with the constraint 



'- 2 + 1 2 )(l + ^-) = ^- + VW. (20) 
a a I a 2 



The dynamics described by Eqs. (|18|) -(|20[) is more complicated in comparison with 
Eqs.(|5|)-(|7]) due to the possibility for the scale factor to have maxima and minima. We 
use the method and notations of jB], where the dynamics of closed universe with the 
scalar field was studied in the absence of the Gauss-Bonnet term. 

Two important regions of the configuration space (a, ip) can be easily found from 
the equations of motion (|18|)-(p0|): the region of possible extrema of the scale factor 
a (for historical reasons it is called the Euclidean one) and the part of the Euclidean 
region where minima of a (i.e. the points of bounce) can appear. 





Figure 1. Regions where points of maximal expansion (gray) and minimal contraction (black) 
are located in configuration space (a, ip). In Fig.f (a) £ = ip 2 , V is asymptotically flat with A < Ao; 
in Fig. 1(b) £ = V 2 , V is asymptotically flat with A > Ao; in Fig. 1(c) £ = ip 2 , V — m 2 ip 2 /2; in 
Fig. 1(d) £ = ip 4 , V = m 2 tp 2 /2. All plots are symmetrical with respect to transition ip — > —-0. 



Eq.(pOl) shows that the Gauss-Bonnet term does not change the configuration of 
the Euclidean region. From the other hand, curves separating the points of maximal 
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expansion (d = 0, a < 0) and the points of minimal contraction (d = 0, a > 0) change 
significantly. Substituting the condition d = into and expressing ip from (p0|) we 
obtain the following equation for the separating curve 

_ 1 + VM + «W_ J< + + = . (2!) 

We start with the case £ = f/> 2 . As for the flat metrics, this case is exceptional 
one. Zero modulus potential leads to asymptotic (2) with constant scale factor 
a = cio = v— 125. The separating curve is simply the line a = do- All the terms 



in constraint (|20|) tend to a constant value on this asymptotic and, therefore, only an 
asymptotically flat potentials could not destroy it. Small enough A corresponds to the 
situation of Fig. 1(a), the nonsingular regime is preserved. This picture changes if A 
exceeds the value A = (2 + \/3)/(4|5|) (Fig. 1(b)). In this case the regime a — > a 
disappears and a trajectory which describes a contracting universe can either fall into a 
singularity or pass trough a bounce and enter an expanding phase. Potentials unbounded 
from above always violate the nonsingular asymptotic (2). 

Now we consider the important case V = m 2 ip 2 /2 in a more detail (Fig. 1(c)). 



The equations (18[)-(20[) were solved numerically using the method of integration 
over an additional parameter developed in Ref. ||. The maximal expansion point was 
chosen as an initial data for the Cauchi problem, unknown values were obtained from 
the constraint equation (p0|). 

If we start from the point with d = 0, initial data corresponding to solutions 
experiencing at least one bounce forms the system of quasiparallel zones in the initial 
condition space (a, ip). This system is rather similar to those studied recently in the 



closed FRW models without the Gauss-Bonnet term ||13|| , but has several important new 
features. 

First of all, it is clear from the separating curve configuration that if the initial 
scale factor is less than \J — 125, the trajectory starts from the bounce itself. So, the 



1-st bounce interval is adjacent to the axis a = in contrast to the results of Ref. [13 
Another important feature can be seen from the Fig.2(a). We plotted cross- sect ions 
ip = of the bounce intervals depending on the mass m. If m < 0.6 we can see the 
full set of intervals: trajectories from 1-st interval have a bounce with no ?/>-turns before 
it, trajectories which have initial point of maximal expansion between 1-st and 2-nd 
intervals fall into a singularity after one ?/>-turn, those from 2-nd interval have a bounce 
after 1 ?/>-turn and so on. If m becomes bigger than 0.6, the 2-nd interval contains 
trajectories with 2 -0-turns before bounce, the space between 1-st interval (which is now 
the product of two merged intervals) and the 2-nd one contains trajectories falling into 
a singularity after two ^-turns. There are no trajectories going to a singularity with 
exactly one ^-turn. 
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Figure 2. The ip = cross-sections of zones in initial condition space (a,ip) leading to bounce 
for V = to 2 V> 2 /2. In Fig.2(a) £ = V> 2 , in Fig2(b) £ = cosh(^) 



This process of interval merging continues with the modulus field mass m increasing. 
The picture is similar to the case without the Gauss-Bonnet term but with gently sloping 
potential ||15|| ) where the process of interval merging also takes place. In both cases, the 
interval formed by merging contains very chaotic trajectories which can not be described 



in a way similar to [13]. A significant part of trajectories starting from such an interval 
may oscillate during very long time without falling into a singularity or leaving the 
interval. In Fig.3. we present the number of trajectories having at least 50 oscillations 
of modulus field versus the modulus field mass. The initial scale factor varies with step 
0.02 in the range of scale factors from 1 to 5 units (total number of trajectories for each 
value of modulus mass is equal to 200). We ignore trajectories which describe a long- 
time expansion (the calculations stopped when running scale factor was 10 times greater 
than the initial one), so, all counted trajectories represented the chaotic oscillations. A 
special effort was done [|T7]] to construct very chaotic solutions in a closed FRW model 
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without the Gauss-Bonnet term, and trajectories having at most 12 oscillations of scalar 
field were found. We choose a much bigger value (50 oscillations), so the existence of 
such trajectories tells us that a qualitatively more complicated chaos exists in the system 
(|I8|)-(|2"0|). If m < 0.7 the measure of such trajectories is so small, that they do not appear 
on our grid. For bigger m this measure grows rather rapidly. For the mass exceeding the 
first merging value 1.3 times more than 30% of trajectories having the point of maximal 
expansion in ip = cross-section of the first interval are very chaotic. So, the regime of 
chaotic oscillation for modulus field masses exceeding the first merging value is much 
more significant than for the case without the Gauss-Bonnet term. 

This picture does not change qualitatively if we switch on the self-interaction of the 
modulus field. For the pure self-interacting potential V(i/)) = the corresponding 
critical value of A leading to the first bounce interval merging is equal to ~ 1.2. 



80 
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Figure 3. The number of trajectories avoiding a singularity during first 50 oscillations of the 
modulus field. Total number of trajectories is equal to 200, initial scale factors are located in the 
range of two first bounce intervals 



When £(V0 grows faster than ijj 2 , there are no solution with asymptotically constant 
scale factor even with the vanishing modulus field potential. For £ ~ ip n+ the 
asymptotic solution is 



/ \ m + 3. ., n-i , ,. i n-i 

a {t) — \ v 1 + ( — 4nd) «+i t ™+i 

V 2n 
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Tp(t) = \-^inr-) t~ i ( 22 ) 



for £ ~ expd^l) the solution was found by Easter and Maeda in ||. In both cases it 
describes the universe with growing scalar factor and modulus field. The Gauss-Bonnet 
term tends to zero and, therefore, at some point it becomes unimportant in comparison 
with arbitrary non-vanishing modulus field potential. After that, the evolution of the 
universe obeys the well-known rules for FRW dynamics with the massive scalar field 
which inevitably lead to recollapse. 

Now it is time to outline the bouncing properties in this case. The separating curve 
looks like plotted in Fig. 1(d). But the if) = cross-section of bouncing intervals is not 
changed much in comparison with the £ ~ ip 2 case. Again, there is the first interval in 
the range of very small scale factors, which is now becoming thinner with the growing 
modulus field. The merging of the bouncing intervals leading to significant chaotisation 
of trajectories is also takes place (Fig.2(b)). 

Additionally we conclude from Fig. 2(b) that though the asymptotic (|2~2"D for zero 
modulus potential is regular, taking into account the classical evolution of the Universe 
trough the maximal expansion point back to the large curvature regime makes the 
dynamics chaotic with the strong dependence on the initial condition. Similar to the 
£ = ip 2 case we have no regular nonsingular solutions, but there are infinite number of 
unstable chaotic trajectories escaping a singularity for an arbitrary long time. 



4. Conclusions 



In this paper we investigated the cosmological behavior of the model with the second 
order curvature corrections based on the action (§) for flat and closed FRW metrics. 
A particular attention was paid to the possibility of nonsingular past asymptotics. We 
confirm recent claims flQj , [TT| , |P2]1 that such regimes exist at least for functions with 
£" > and present corresponding formulas for power-law and exponentially steep £(V0- 
The previously known case £ ~ ip 2 appears to be in some sense the exceptional one. It 
corresponds to DeSitter flat universe or closed universe with the constant scale factor 
(see Eqs.(2),(H)). Steeper £(-0) lead to power-low dependence of the Hubble parameter 
upon time for the flat metric and to bounce nonsingular solutions for the closed one (a 
particular example with exponentially steep £(V0 for the latter case was also studied by 
Easther and Maeda in The important point also mentioned in above cited papers 
is that these asymptotic regimes do not require a fine-tuning of initial conditions. 

This optimistic view on the nonsingular regimes changes, however, if we take 
into account possible potential term for modulus field, which naturally arises in the 
description of SUSY breaking. We conclude that a modulus potential destroys the 
nonsingular asymptotics in most cases we have studied. For a flat universe the only case 
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with £" > which allows nonsingular regime and a massive modulus field is £ ~ ip 2 . 
Even in this case a steeper modulus potential kills the nonsingular asymptotic. For 
£ ~ ip 4 it can survive only for asymptotically flat potentials. Steeper £ can not prevent 
from falling into a singularity with an arbitrary modulus field potential not tending 
to zero for — > oo. Our numerical simulations show that when the nonsingular 
regime described above becomes impossible, all the cosmological trajectories go towards 
a singularity and no new nonsingular regime appear. 

For a closed universe the possibility of regular nonsingular asymptotic is even more 
restricted. The above mentioned asymptotic can survive only in the exceptional case 
£ ~ ip 2 and an asymptotically flat potential. But now a trajectory which does not 
belong to this regime may either fall into a singularity or experience bounce, depending 
on the initial conditions. At the expanding phase after bounce the Gauss-Bonnet term 
become irrelevant and we can use the results obtained for the scalar field in Einstein 
theory. Since for the closed metric all the trajectories must have their point of maximal 
expansion, a bouncing universe will begin to contract at some point and, at the second 
time, fall into a singularity or have a second bounce etc. So, we return to the chaotic 
picture similar to one described for the massive scalar field coupled with the Einstein 
gravity. This chaotic regime is very sensitive to the initial condition and does not contain 
any stable nonsingular trajectory. 

For steeper the regime without a modulus potential is a bounce itself, so 

arbitrary potential at the after-bounce stage will dominate the Gauss-Bonnet term and 
lead to recollapse and chaotisation of trajectories. 

Since the particular form of second order curvature corrections arising from the 
string theory include exponentially steep with coming from SUSY breaking 

exponentially steep modulus potential, nonsingular solutions presenting in the string 
gravity with V(ip) = disappear when modulus potential is taken into account. 
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